Quantum-Limited Mirror-Motion Estimation 
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We experimentally demonstrate optomechanical motion and force measurements near the quan- 
tum precision limits set by the quantum Cramer- Rao bounds (QCRBs). Optical beams in coherent 
and phase-squeezed states are used to measure the motion of a mirror under an external stochastic 
force. Utilizing optical phase tracking and quantum smoothing techniques, we achieve position, 
momentum, and force estimation accuracies close to the QCRBs with the coherent state, while esti- 
mation using squeezed states shows clear quantum enhancements beyond the coherent-state bounds. 
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The advance of science and technology demands in- 
creasingly precise measurements of physical quantities. 
The probabilistic nature of quantum mechanics rep- 
resents a fundamental roadblock. Over the last few 
decades, the issue of quantum limits to precision mea- 
surements has been a key driver in the development 
of quantum measurement theory |2[. With the re- 
cent technological advances in quantum optical, electri- 
cal, atomic, and mechanical systems, quantum limits 
are now becoming relevant to many metrological appli- 
cations, such as gravitational- wave detection [3], force 
sensing mametometry [5], clocks [6], and biological 
measurements [7|. 

It is now recognized that quantum detection and es- 
timation theory [8| provides the appropriate framework 
for the definition and proof of quantum measurement lim- 
its. For parameter estimation and the mean-square er- 
ror (MSE) criterion, a widely studied quantum limit is 
the quantum Cramer- Rao bound (QCRB) [1, [9|. For 
gravitational-wave astronomy and many other sensing 
applications, the estimation of time- varying parameters, 
commonly called waveforms in the engineering literature, 
is more relevant. QCRBs for waveform estimation were 
recently derived in Refs. 0, although there has not yet 
been any comparison of the waveform QCRBs with ex- 
perimental results to demonstrate their relevance to cur- 
rent technology. 

Quantum estimation of an optical phase waveform was 
recently demonstrated experimentally 11, [l2| using an 



optical phase tracking method that measures t he p hase 
via homodyne detection with feedback control [13|, fol- 
lowed by smoothing of the data Q. These experi- 
ments demonstrate improvements over heterodyne mea- 
surements, causal filtering [ll|, and coherent-state opti- 
cal beams when squeezed light is used [12|, but no com- 
parison with the QCRBs was made to test the optimality 
of the experimental techniques. 



In this Letter, we report an experiment that applies 
the tracking and smoothing techniques to optomechani- 
cal motion sensing. We use optical probe beams in coher- 
ent and phase-squeezed states to measure the motion of a 
mirror under an external stochastic force and then com- 
pare the smoothing errors with the waveform QCRBs. 
This is the first time to our knowledge that experimental 
results have been compared with the waveform QCRBs. 
Through the comparison, we are able to demonstrate the 
near-optimality of our measurement method in the case 
of coherent states. The squeezed-state results are further 
away from the QCRBs but still show clear enhancements 
over the coherent-state bounds. Despite our focus here 
on a classical mechanical system, our methods can also be 
applied to purely quantum systems 1Q|, [H | , making our 
methods potentially useful for a wide range of quantum 
sensing applications [l-Q- 
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FIG. 1. (Color online) (a) Schematic of mirror- motion es- 
timation, (b) Experimental setup. LO: Local Oscillator, 
RF: Radio Frequency, Ti:S: Titanium Sapphire laser, AOM: 
Acousto-Optic Modulator, EOM: Electro-Optic Modulator, 
SHG: Second Harmonic Generator, OPO: Optical Parametric 
Oscillator, FPGA: Field-Programmable Gate Array. 

Figure [TJa) shows a schematic of our experiment. 
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where the mirror motion is approximated as a mass- 
spring-damper system. The mirror, driven by a stochas- 
tic force, is illuminated by a probe beam in a coherent 
state or a phase-squeezed state. The motion of the mir- 
ror shifts the phase of the probe beam. We measure this 
phase shift adaptively by homodyne detection (optical 
phase tracking) |llMl3|. and estimate the mirror motion 



14. 



from the optical phase measurements 

Optical phase tracking allows us to linearize the mea- 
surement results y{t) as 



y{t) = ^{t)^z{t), 



(1) 



where ip{t) is the optical phase shift and z{t) is a noise 
term depending on the optical beam statistics |llUl2Ul5| . 
The phase shift (p{t) of the probe beam is caused by the 
mirror position shift q{t) as 



(p{t) = {2ko cos 0)q{t), 



(2) 



where ko cos is the wave- vector component parallel to 
the mirror motion and is the reflecting angle as shown in 
Fig. [11(a), fixed at 7r/4. We estimate the mirror position 
q{t)^ momentum p(t), and external force f{t) from the 
measurement results y{t). q{t)^ p{t)^ fif)i ^"^^ ^(^) 
assumed to be zero-mean stationary processes. 

Under the linear approximation, the optimal estimate 
of the mirror position is a weighted sum of the measure- 
ment results given by q'{t) = drJq{t — r)y{r)^ where 
Jq{t) is a linear filter and prime indicates an estimate. 
Estimates of momentum p\t) and external force f\t) 
are similarly defined. The integration limits are approxi- 
mated as ±00 because we use data long before and after 
t to obtain the estimates at the intermediate time t via 
smoothing The optimal position filter Jq{t) is ob- 

tained by minimizing the MSE 11^ = {[q\t) - q{t)]'^) , 
which is averaged over the probability measures for z{t) 
and q{t) {Hp and Hf are similarly defined). The optimal 
filters and the minimum MSEs are calculated by moving 
to the frequency domain il5il . The minimum MSEs H^^" 
{x = q^p^ f) are given by |14l4l6| 
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where Sxioo) {x = q^p^f^z) is a spectral density defined 
as Sx{i^) '= dr {x(t)x{t + r)) e*^'^, g^pxi^) is a trans- 
fer function that relates the optical phase shift (f to the 
target variables {x = /) by (f{(jo) = gcpx{^)x{uj)^ with 
the tilde indicating a Fourier transform. 

We now consider the QCRBs on the MSEs. The wave- 
form QCRBs are derived from the quantum properties of 
the probe beams and prior statistics of the target system 
(mirror motion) and do not depend on the measurement 
and post-£rocessing method. The QCRBs for our situa- 
tion are llQl 
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where Sai{^) is the spectral density of the probe-beam 
photon flux. Comparing Eq. (|3]) with Eq. (j4]), we find 
that 4.Sai{oj) = l/Sz{oj) is required for to match 

the QCRBs. This means that, to attain the QCRBs, 

(i) the probe beam should be in a minimum-uncertainty 
state with respect to the phase and the photon flux, and 

(ii) the measurement noise z{t) should consist of intrinsic 
phase noise only. 

Our experiment uses broadband phase-squeezed states, 
including coherent states as the small-squeezing limit. 
The noise term z{t) in the normalized homodyne outputs 
can be written in a quadratic approximation [lH, [l5| as 

(5) 
(6) 
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where rm(rp) is the squeezing (anti-squeezing) parameter 
(^p > > 0), a is the coherent amplitude of the probe 
beam, is the steady-state MSE of the optical phase 
estimate in the real-time feedback l oop (a^ <C 1). ^sq is 
called the effective squeezing factor [125, which takes into 
account the anti-squeezed amplitude quadrature as well 
as the squeezed phase quadrature. The noise spectral 
density Sz{uj) and the photon- flux spectal density Sai{^) 
are 
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Here we assume that the bandwidth of squeezing is broad 
compared to the bandwidth of system parameters, but 
not too large so that the photon-flux fluctuations do not 
diverge (see Supplemental Material [HI). 

The necessary condition to reach the QCRBs is now 
given by e^^p = 1/^sq- For coherent states (rm = 
Tp = and ^sq = 1), this condition is always sat- 
isfied, so QCRB-limited estimation is possible within 
the quadratic approximation. On the other hand, 
the squeezed-state QCRB is attainable only if (i) the 



squeezed state is pure (e^ 



e^^"^) and (ii) the opti- 



cal phase tracking works well enough such that 0. 
Thus, in a real experimental situation, the squeezed-state 
QCRB is more difficult to reach than the coherent-state 
QCRB. We emphasize however that our estimation re- 
sults are still comparable to the squeeze-state QCRBs 
and better than the coherent-state bounds. 

Figure [Hb) shows our experimental setup. A 
continuous-wave Titanium Sapphire laser is used as a 
light source at 860 nm. Phase-squeezed states aregener- 
ated by an optical parametric oscillator (OPO) [12I, Il7|. 
The OPO is driven below threshold by a 430 nm pump 
beam. Optical sidebands at ± 5 MHz are used as a car- 
rier beam generated by acousto-optic modulators |llUl2|. 
To avoid experimental complexities, the pump power is 
fixed at 80 mW, producing squeezing and anti-squeezing 
levels of -3.62±0.26 dB and 6.00±0.15 dB. The effective 
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squeezing factor, ^gq, varies from —3.28 dB to —3.48 dB 
depending on the probe amplitude. To make a coherent 
state, we simply block the pump beam. 

A mirror (12.7 mm in diameter, 1.5 mm in thickness, 
0.444 g in weight) is attached to a piezoelectric trans- 
ducer (PZT, weighing 0.432 g). We assume the mass of 
this PZT-mounted mirror to be m = (0.444 + 0.432/3) g 
= 5.88 X 10-^ kg from the uniformity of the PZT pj- 
The transfer function of the PZT-mounted mirror (the re- 
lation of applied voltage to actual position shift) is mea- 
sured before the estimation experiments. We use this 
transfer function to construct optimal filters and calcu- 
late the QCRBs [lij. 

In the estimation experiments, the PZT-mounted mir- 
ror is driven by an Ornstein-Uhlenbeck process. This sig- 
nal is generated by a random signal generator followed by 
a low-pass filter with a cutoff frequency of A = 5.84 x 10^ 
rad/s. We drive the PZT within the linear response range 
so that the external force f{t) is proportional to the sig- 
nal. Thus the external force f{t) is also an Ornstein- 
Uhlenbeck process given by 



dfjt) 
dt 



-\f{t)+w{t), 



(8) 



where w{t) is a zero- mean white Gaussian noise satisfying 
{w{t)w{r)) = i<i6{t — r). In the experiment, we set hz = 
1.67 X 10^ s-^ 

A fraction of the laser beam is used as a local oscillator 
beam, which is optically mixed with the probe beam at 
a 1:1 beam splitter for homodyne detection. The over- 
all efficiency of the detection is 87% jl5|. The homo- 
dyne output is demodulated and recorded with an os- 
cilloscope. The measured data are post-processed us- 
ing a computer to produce the estimates. The demodu- 
lated homodyne output is also processed by a field pro- 
grammable gate array (FPGA) for the real-time feedback 
based on Kalman filtering, which approximates the mir- 
ror motion as a mass-spring-damper system [l^. Note 
that we use this approximate model only for the real-time 
feedback, not for the estimation. In the experiment, we 
have another low-gain, low-frequency feedback loop to 
prevent environmental phase drift. 

Figure [2] shows one of the time-domain results for 
the mirror-motion estimation with phase-squeezed states. 
The black lines are the signals to be estimated (for the 
evaluation, see Supplemental Material (lH). The exter- 
nal force / is an Ornstein-Uhlenbeck process given by 
Eq. ([8]). The periodic oscillations of q and p arise from 
the mechanical resonance of the PZT-mounted mirror, 
the frequency of which is 1.76 x 10^ rad/s [l5|. The red 
lines are the estimates, which agree well with the sig- 
nals. This 1 ms long data are obtained with a sampling 
frequency of 10 MHz, and are repeated 300 times to eval- 
uate the MSEs. 

We perform mirror-motion estimation with probe 
beams in the coherent state and the phase-squeezed state. 
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FIG. 2. (Color online) Time-domain results for (q) posi- 
tion, (p) momentum, and (/) external force, respectively, with 
|q;|^ = 6.24 x 10^ s~^ and the probe beam in a phase-squeezed 
state. The black lines are the signals to be estimated. The 
red lines (gray lines in print) are the estimates. 



each with four different amplitudes. Figure [3] shows the 
|ap dependence of the MSEs of the position, momen- 
tum, and external force estimation. Figure [3] shows three 
key results. First key result: Experimental results agree 
well with the theoretical predictions (traces i and iii). 
The small discrepancies may be attributed to the low- 
frequency noise due to environmental phase drift, and 
slight changes of the mirror properties (e.g., the reso- 
nant frequency) during the experiment. Second key re- 
sult: The experimental results are close to the wave- 
form QCRBs. In particular, the experimental results 
for coherent states (green circles) are very close to the 
coherent-state QCRBs (traces ii). The closeness (i.e., 
relative differences between the experimental MSEs and 
the coherent- state QCRBs) is quantified as 28 ± 12%, 
15 ± 6%, and 11 ± 6% on average for the position, mo- 
mentum and force estimates, respectively. The small dif- 
ferences between the prediction curves (traces i) and the 
coherent-state QCRBs (traces ii) are attributed to the 
imperfect detection efficiency. The experimental results 
of squeezed states (red diamonds) are also comparable to 
the squeezed-state QCRBs (traces iv), although the gaps 
are larger due to the impurity of the squeezed states. 
Third key result: The experimental results for squeezed 
states show clear quantum enhancement, mostly over- 
coming the coherent-state QCRBs. The quantum en- 
hancements (i.e., relative reduction of MSEs compared 
to the coherent-state QCRBs) are quantified as 15 ± 8% 
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FIG. 3. (Color online) Experimental and theoretical MSEs 
of the (q) position, (p) momentum, and (/) external force, 
plotted versus the probe amplitude squared, The green 

circles are the results for coherent states, and the red dia- 
monds are those for phase-squeezed states. The green solid 
curves (traces i) are simulated prediction curves of the esti- 
mates, which were calculated by considering the experimental 
imperfections. The green dot-dashed curves (traces ii) are the 
coherent-state QCRBs. The red solid lines (traces iii) are the 
simulated prediction curves for a phase-squeezed probe beam, 
where we use the quadratic approximation as in Ref. p^ . 
The red dot-dashed curves (trace iv) are the squeezed-state 
QCRBs. 



and 12 ± 2% on average for the position and momentum 
estimates, respectively. The force estimate at the highest 
probe amplitude is slightly worse than the coherent-state 
QCRB, which should be due to the low- frequency noise 
from the environment. Note that we still observe quan- 
tum enhancement of the force estimation (except the esti- 
mate at the highest probe amplitude), which is quantified 
as 12 ± 2% on average. 

In conclusion, we have experimentally demonstrated 
quantum-limited mirror-motion estimation via optical 
phase tracking. Our experiment reveals that the 
coherent-state QCRB is almost attainable by our ex- 



perimental method. Although the squeezed-state QCRB 
turns out to be more difficult to reach because of the im- 
purity of the squeezed states, quantum enhancement be- 
yond the coherent-state QCRB is clearly observed. These 
results demonstrate the potential of our theoretical and 
experimental methods for future quantum metrological 
appHcations. 
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Supplemental Material for 
Quantum-Limited Mirror-Motion Estimation 

EXPERIMENTAL DETAILS 

In this section, we will describe the experimental de- 
tails. Figure m shows our experimental setup [12|. A 



RF source 

— >\o mixer 



Signal 



continuous-wave Titanium Sapphire laser was used as a 
light source at 860 nm. Phase squeezed states were gen- 
erated by an optical parametric oscillator (OPO) of a 
bow-tie shaped configuration with a periodically polled 
KTiOP04 crystal as a nonhnear optical medium (l7| . 
The OPO was driven below threshold by a 430 nm pump 
beam, generated by another bow-tie shaped cavity that 
contains a KNbOs crystal. The free spectral range and 
the half width at half maximum of the OPO were 1 GHz 
and 13 MHz respectively. Optical sidebands at ± 5 MHz 
were used as a carrier beam generated with acousto-optic 
modulators 



III, ll2|. 



Note that these optical sidebands 
are within the OPO's bandwidth. To avoid experimental 
complexities, the pump power was fixed to 80 mW giv- 
ing squeezing and anti-squeezing levels of —3.62 zb 0.26 
dB and 6.00±0.15 dB respectively. The effective squeez- 
ing factor, ^sq, varied from —3.28 dB to —3.48 dB ac- 
cording to the probe amplitude. Note that takes 
into account of the anti-squeezing quadratures mixing 
in the measurement, which cannot be neglected for rel- 
atively high squeezing levels. It is a trade-off between 
enhancement from the squeezed quadratures and degra- 
dation from the anti-squeezed quadratures, revealing an 
optimal squeezing level [l2| . The optimal squeezing level 
differs for each amplitude but the difference is mi- 
nor for our experimental conditions. Since the gener- 
ated phase squeezed state becomes less robust for higher 
pumping levels due to the complex locking system, we 
chose a slightly lower pumping level and did not change 
it for each \a\. For comparison to phase squeezed states, 
we also used coherent states as a probe by simply block- 
ing the pump beam. 

The mirror mounted on a piezoelectric transducer 
(PZT) was driven by a signal that follows the Ornstein- 
Uhlenbeck process. This signal was generated with a ran- 
dom signal generator followed by a low-pass filter with a 
cutoff frequency of A = 5.84 x 10^ rad/s. 

A fraction of the laser beam was used as a local oscilla- 
tor (LO) beam which was passed through a spatial-mode 
cleaning cavity (not shown in Fig. |4]) to increase mode 
matching with the probe beam. The probe beam and 
the LO beam are optically mixed with 1:1 beam splitter 
for homodyne detection. The efficiency of the detection 
is shown in Table |T1 The homodyne output was demodu- 
lated and recorded with an oscilloscope for post process- 
ing. 

In the feedback loop, the LO phase is modulated ac- 
cording to the estimated phase. The modulation was 




FIG. 4. Experimental setup. Ti:S: Titanium Sapphire laser, 
LO: Local Oscillator, RF: Radio Frequency, AOM: Acousto- 
Optic Modulator, EOM: Electro-Optic Modulator, SHG: Sec- 
ond Harmonic Generator, OPO: Optical Parametric Oscilla- 
tor, FPGA: Field Programmable Gate Array. 



TABLE L Efficiency of the detection. 



Photo diode quantum efficiency 


0.99 




Interference efficiency (Visibility) 


0.965 


(0.982) 


Propagation efficiency 


0.981 




Electrical circuit efficiency (Clearance) 


0.924 


(11.2 dB) 


Overall efficiency 


0.871 





performed with a waveguide type electro-optic modula- 
tor (EOM). The real-time phase estimate used for feed- 
back was processed with a field programmable gate array 
(FPGA). The delay of our implemented feedback filter 
was around 400 ns, which is small enough for our cur- 
rent experimental parameters. Note that we have an- 
other low-gain, low-frequency feedback loop to prevent 
environmental phase drifting. 



MODELING THE MIRROR MOTION 

In this section, we will explain modeling the mirror 
motion. First, we will consider how to evaluate mass of 
a PZT-mounted mirror. Then, we will explain transfer 
function of the PZT-mounted mirror, and the evaluation 
of true signals to be estimated. Finally we will describe 
the mirror motion functions. 



Mass of a mirror attached to a PZT 

In our experiment, a multilayer PZT (AE0203D04F, 
NEC/Tokin) of 3.5 mmx4.5 mmx5.0 mm in size weigh- 
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ing 0.432 g was used. A mirror, 12.7 mm in diameter, 
1.5 mm in thickness, weighing 0.444 g was attached to 
the PZT with an epoxy-based adhesive. The mass of the 
mirror attached to the PZT was evaluated as fohows. 

Let the mass of the PZT and mirror be Mp and Mm, 
respectively. Assume that the mass of the PZT is uni- 
form, and that the displacement is proportional at all 
points. 



Al = — AL. 



(9) 



Here, the original length of the PZT is Lq, the overall 
displacement is AL, and the displacement at point / (0 < 
/ < Lo) is AL Then, the kinetic energy may be calculated 

as 



M„ 



1 2 



dt 



Hence, we assume that m = 
0.432/3) g = 5.88 x lO""^ kg. 



(AL) 



(10) 



Mp/3 = (0.444 




OU cutoff 



50 100 
Frequency (kHz) 
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FIG. 5. Transfer functions of the PZT-mounted mirror, gain 
(a) and phase (b) . Black solid lines show the measured trans- 
fer function Tq(uj) referred to as the detailed model. Red 
dashed lines show the fitted transfer function of the mass- 
spring-damper system Tq^'^^uj) referred to as the nominal 
model. The green dot-dashed line shows the cutoff frequency, 
A/27r, of the Ornstein-Uhlenbeck (OU) process used in the 
experiment. 



Transfer function of the PZT-mounted mirror 

Next, we will focus on modeling the transfer function 
of the PZT-mounted mirror. The mass-spring-damper 
model is referred to as the nominal models which is a 
simplified model that describes the essence of the tar- 
geted system. On the other hand, a model which best 
describes the targeted system is referred to as the de- 
tailed model. The detailed model would be the closest 
measurable model of the targeted system. We used this 
detailed model to construct optimal filters and calculate 
the QCRBs, while we used the nominal model to realize 
real-time feedback control. 

We used a Mach-Zehnder interferometer and a net- 
work analyzer to measure the transfer function of the 
PZT-mounted mirror, To(a;), referred to as the detailed 
model. The black solid lines in Fig. [5] show the mea- 
sured results. The red dashed lines in Fig. [5] show the 
fitted transfer function of the nominal model T(5^°"^(cj) ex 
1/ (— mcj^ + imcJ7 + ml^^) where 7 is the damping co- 
efficient and is the mechanical resonant frequency. 
The fitted parameters were Vt = 1.76 x 10^ rad/s and 
7 = 7.66 X 10^ rad/s. 

Note that the external force driving the mirror is gen- 
erated according to the Ornstein-Uhlenbeck process. The 
cutoff frequency of this process was set to A = 5.84 x 10^ 
rad/s, which is indicated as a green dot-dashed line in 
Fig. [5l The nominal model is good enough to construct 
the real-time feedback filter for the experimental condi- 
tions. 



Evaluation of true signals 

In order to evaluate estimation errors, we need to know 
the true position, momentum and external force that are 
to be estimated (referred to as the target position, target 
momentum, and target force). We use the full range of 
the detailed model Tq{uo) to calculate these target posi- 
tion momentum p, and external force /. In the mir- 
ror motion estimation experiment, we record the volt- 
age V{t) that drives the PZT-mounted mirror. From 
V(t), To(co'), and the sensitivity of the photo detector 
G = 6.96 X 10^ V/m, we calculate the target position as 



q{t) =T- 



G 



(11) 



where T {T ^) denotes the (inverse) Fourier transform. 
We use this result to calculate the target momentum. 



(12) 



The voltage applied to the PZT-mounted mirror is within 
the linear response range so that the target force may be 
calculated as 



j{t) = pV{t), 
where /3 = 2.04 x lO'^ N/V. 



(13) 
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Mirror motion functions 

Mirror motion functions gij{oo) {i^j = q^p^f^^) are 
defined such as q{uj) = gqf{uj)f{uj) {i = q and j = /) 
where a tilde indicates the Fourier transform. The mir- 
ror motion functions are necessary to derive the optimal 
filters and the QCRBs. Note that the definition leads to 
9ij{^)9jk{^) = 9ik{^) and [gij{uj)]~^ = 9ji{^)' 

From Eqs. (pT|) and (fT3]) . the function Qqfioo) is given 
as, 



G(3 ' 



(14) 



As denoted in the main text, the phase shift of the probe 
beam is proportional to the position shift as (p{t) = 
{2ko cos 0)q{t). Then, the other relevant mirror motion 
functions are derived as follows: 



9^q{^) =2/^0 cos6>, 

9^p{^) =9^q{^)9qp{^) 







(15) 


2ko cos/ 
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(16) 




imuj 




2ko cost 


G(3 ' 


(17) 



9^f{^) =9m^^)9qf{^) = 2A;ocos 
where we use gqp{oo) = l/{imuj). 



OPTIMAL LINEAR FILTER AND LEAST MEAN 
SQUARE ERROR 

In this section, we derive the optimal linear filters 
which minimize mean square errors (MSEs) [l6| . We will 
explain the position estimate q^t) and the least position 
MSE n^^" as an example. The estimates and MSEs for 
momentum and force can be derived similarly. 

First, let's consider the normalized output of the ho- 
modyne detection 12, 13| . 



rj{t) = sinMt) - ^'{t)] + 1^ V'^), (18) 



Rsq{t)=sm^[ip{t)-^'{t)]e^-- 



-2r„ 



(19) 



Here rm (^p) is the squeezing (anti-squeezing) parameter 
(^p > > 0), |a| is the coherent amplitude of the probe 
beam, v{t) denotes white Gaussian noise with a flat spec- 
tral density of 1, and (p'{t) is a real-time phase estimate 
used for the feedback control. This homodyne output 
can also be applied to coherent states by simply putting 
Rsq = 1. Following the quadratic approximation shown 
in Ref. [12| gives a good approximation of the homodyne 
output as 



Here, z{t) is a white Gaussian noise as, 

{z{t)) = 0, 



{z{t)z{r)) 
R. 



sq 



4|a|2 

^2 ^2r, 



S{t-T), 



«q - a^^- + (1 - (T^)e- 



-2rn 



(21) 

(22) 

(23) 
(24) 



^sq is called the effective squeezing factor [12], which 
takes into account the anti-squeezed amplitude quadra- 
ture as well as the squeezed phase quadrature. 

By adding the real-time phase estimate (p'{t) (which is 
measured in the experiment as well as r]{t)) to r]{t)^ we 
obtain the (modified) measurement result y{t)^ 



y{t) = 7^{t) + ^'{t)c^^{t) + z{t). 



(25) 



The linear estimate of position, q\t), is given as a 
weighted sum of this y{t)^ 



+ 00 



drJqit- T)y{T), 



(26) 



where Jq(t) is a linear position filter. Fourier transform 
of the estimate is calculated as, 

q'{uj) =Jq{uj)y{uj) = Jq{uj) [(f{uj) + z{uj)] 

= Jq{^) [9^q{^)q{^) + ^{UJ)] . (27) 

We define a two-time covariance Sg(t, t + r), 

S,(t, t + T):= {[q'{t) - q{t)] [q'{t + r) - q{t + r)]) . 

(28) 

Note that we stick to steady-state so that T,q{t^t + r) is 
determined by only r. The Fourier transform of Eg(t, t + 
r) is defined as. 



/ + 00 
(iTS,(t,t + r)e'"". 
-oo 

MSE of the position estimation, 11^, is given as as, 



n, := 



J — C 



27r 



(29) 



(30) 



Our aim is to derive the filter Jq{oj) minimizing 11^ and 



obtain the least Hq. 



Let's focus on Cq{uj) because Uq is minimized by min- 
imizing Cq{uj) at all the uj. After some algebra, we find 
the following: 



Cq{uj) = Jq{(jj)g^q{(jj) - 1 Sq{(jj) + Jq{uj) S 



(31) 



where Sk{oj) is a spectral density defined as Sk{oj) 



(20) J+^^^^^(^)^(^^^)^g^^r 



q^z). By setting 
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dCq{(jo) / dJq{(jo) = 0, we obtain the optimal position filter 



g*{u)Sq{uj) 



Accordingly, the least MSE is derived as, 



TTinin 



2^ \Sq{uj) S,{uj) 



(32) 



(33) 



The other optimal filters and MSEs for p and / can be 
obtained by changing the subscript q to p or f. 

The spectral densities Sk{oj) {k = q^p^f^z) in our ex- 
periment are obtained as follows: First, Sz{oj) is easily 
obtained from Eq. (|22]) , 



'sq 



4|a 



12 • 



(34) 



The external force f{t) obeys the Ornstein-Uhlenbeck 
process. 



df{t) 
dt 



{w{t)w{r)) n5{t-T). 
Thus, Sf{uj) is given as. 



(35) 
(36) 

(37) 



Other spectral densities can be calculated by using the 
relation Si{uj) = \gij{uj)\'^Sj{uj). From Eq. (fT4|) we ob- 
tain. 



G/3 



w2 + A2 ' 



(38) 



Sp{oj) = \gpf{uj)\'^Sf{uj) = \gpq{uj)gqf{oj)\'^Sf{uj) 



G/3 



(39) 



PHOTON FLUX FLUCTUATION 

In this section, we will derive the spectral density of the 
photon flux fluctuation S'a/(w) and discuss the validity 
of the approximation used in the main text, S'a/(w) « 

In order to calculate the photon flux fluctuation, we 
use an annihilation operator for an electromagnetic fleld, 
a{t), which satisfies the commutation relation jl8l |. 

[a{t),a^{t')]=S{t-t'), (40) 
[a{t),a{t')]= [a\t),a\t')]=0. (41) 



Photon flux I{t) and the mean photon flux Iq are given 
as, 



I(t) =a\t)a{t), 



(42) 
(43) 



We define the Fourier transform of an annihilation oper- 
ator, 



/ + 00 
dta{t) 
-co 



Aujt 



(44) 



Note that this definition leads to d^u) = d{—uj). The 
commutation relation in the frequency domain is derived 
from Eqs. (j40]) and (|4T]) . 

[a(cj), d\u)] = 27r5{uj - u), (45) 
[d{uj),d{uj')] = [a^(cj),a^(cjO] = 0. (46) 
Spectral density of a{t) is given as. 



R{uj) :-- 



dT {a\t)a{t + t)) e''^^ (47) 
^(at(c.i)a(c.)). (48) 



The mean photon flux Iq is obtained by integrating this 
spectral density R{uj)^ 



(49) 



The spectral density of the photon flux fluctuation 
Sai{^) is calculated as. 



+ 00 



dr((/(t)-/o)(/(t + r)-/o))e^' 



duJiduj2duJ3 /-+/ x-f/ x~/ x~/ x\ 
— ^^-^3 — {a^ {-ujija^ {-uJs)a[uj2)ci{oj - 0J3)) 



00 

+ CO 



— i?(w) — / dwidw2-R(wi)-R(w2)- 

27r 27r 

(50) 

To derive 5'a/(^), we have to calculate the fourth order 
moment of an annihilation operator. In our case, how- 
ever, we use a Gaussian state (phase squeezed state), so 
the second order moment will suffice to describe S'a/(^)- 
Let's assume an annihilation operator of the form, 

d{uj) =27i5{uo) \a\ + asq(a;), (51) 
asq(a;) =cia(a;)ai(cj) + cib(a;)a|(-a;) 

+C2a(^)tt2(^) + C2h{^)d\{-Uj), (52) 

where \a\ is a coherent amplitude, d^o^iuo) represents the 
squeezing term ((asq(co')) = 0), di{(jj) and a2(cj) are vac- 
uum modes. Here we set the amplitude as a real value 
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without loss of generality. The expression of Eq. (j52j) is 
valid for any mean-zero Gaussian states including mixed 
states (i.e., squeezed thermal states), as long as the coef- 
ficient Cij{uj) satisfies the following: 



c;juj) =Cij{-uj), 



(53) 



\cMf - \cMf + |c2a(c^)|' - |c2b(c^)|' = 1, (54) 



cL(^)cib(^)+C2a(c^)c2b(c^) 
-Cia(cj)Cib(^)-C2a(^)c2b(^) = 0- 



(55) 



Here these equations are imposed by the property of the 
Fourier transform and the commutation relation (Eqs. 
(pj) and (|i6|) ). 

To describe the photon flux fluctuation of the squeezed 
states, it is useful to define the quadrature operators, 

Xsq(cj) [asq(cj) + ajq(-cj)] , (56) 

Psq(^) [asdi^) - aldi-^)] • (57) 

Since we set the amplitude as a real value, Xgq (Psq) is 
the anti-squeezing (squeezing) quadrature. Photon flux 
spectrum (except the amplitude contribution), squeezing 
spectrum and anti-squeezing spectrum (spectral densities 
of asq(co'), Psq(^) and Xsq{uj)) are given as. 



-oo 27r 

+00 



(alq(wi)asq(w)) , 



(58) 
(59) 
(60) 



Here squeezing and anti-squeezing spectrum satisfy an 
uncertainty principle, R'^^{uj) R~^{uj) > 1/16 {h = 1/2). 
From Eqs. (f5T]) ~ ([60]) . we obtain. 



^sq(^) =7 |cia(c^) ±Cib(c^)r 



C2a(^) ±C2b(^)| , 



Ri^iu;) =i?+(c^) + i?-(c^)--. 
Then, after some algebra, Eq. (|5Q]) is rewritten as. 



(61) 
(62) 



i 



-sqV 

+°° duji 
27r 

27r 



2i?+(c^i+w)i?i('^i) 



2i?-(wi+w)i?-(wi) 



■ (63) 



Next, we will assume that the squeezed state has finite 
bandwidth, and then verify the approximation used in 
the main text. 



Let's consider the squeezing spectrums Rf^{(jo) of the 
standard form 12|, |T3 , 




(Ac^±)' 



(Ac^±)^ 



4i?s-q(0) 



4i?+ (0) - 1 



(64) 



(65) 



where Acj_ (Acj+) is the bandwidth of squeezing 
(anti-squeezing), the second equation ensures that 
R'^J^(jj)R~^{(jj) = 1/16 for all a; when the squeezed state is 
pure. Here we define the averaged squeezing bandwidth 
AcjOq and the squeezing parameters (rm, ^p) at the center 
frequency {uj = 0) as. 



Aujo 



: -(Acj_+Acj+), 

:4i?-(0), 

4<(0). 



(66) 

(67) 
(68) 



In the case of the OPO, Aujq corresponds to the half 
width at half maximum of the OPO. 

By inserting Eq. (f64|) to Eq. (|63]) . we obtain. 



5a/(w) =4|a|2i?+(u;)+J,q 

(gS-P _ l)2(Aw+)3 (1 



02(Aw_)3 



^sq 



(2Au;+)2 



(2Aa;_)2 



27r 



l)Aw+ + (e- 



-2r„ 



l)Au 



(69) 



(70) 



where /gq is the mean photon flux of the squeezing (/q 



\a 



^sq 
^+/sq). 



If the averaged squeezing bandwidth Auq is much 
larger than the system parameters, i.e., Auq ft^X {Q: 
the resonant frequency of the mirror. A: the cutoff fre- 
quency of the external force), we may assume u <C Auj±. 
Note that we implicitly assume that Auq ~ Au- ^ Aa;+, 
which would be justified in our experimental situation 
as described later. The photon flux fluctuation Sai{^^) 
would be approximated to. 



SAi{io)-{\af+^Isci)e 



2r„ 



1 (e^'- 
4 



1)3/2 ^ (1 



1-VT 



-2r„ 



where the parameter ^ ranges from 1 (rp 
1/4 (rp^(^). IfC/sq«|ap, 

S'a/(w) 



(71) 

e-2'--^)3/2 ^ 

(72) 
0) to 



(73) 



Let's consider whether these conditions {Aujq ^ Q^X 
and ^/gq <^ |<^P) are satisfied under our experimen- 
tal situation. The experimental parameters are, e^^p = 
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3.98 (6.00 dB), e-2^-=0.435 (-3.62 dB), ^ = 0.61, 
Aa;+/Aa;_ = 0.435, Q = 1.76 x 10^ rad/s, and A = 
5.84 X 10^ rad/s. The averaged squeezing bandwidth 
AcJo, however, is tricky to determine. As in Ref. [l2| , 
we utiHze only finite bandwidth around a sideband fre- 
quency of 5 MHz. Thus it is not appropriate to define 
the squeezing bandwidth Aujq as an OPO's bandwidth 
(Acjopo = 8.2 X 10^ rad/s ). We should consider an 
effective squeezing bandwidth AcJq^ which is not unnec- 
essarily large, but still satisfies Aujq^ 1], A. 



(In 





' (iv)"|ap =6.24x10^ 






(iii) \af =2.87x10^ / 






fii) \a\' =1.88x10^ / 






(i) lap =1.02x10^ X 




A 


0, Acjqpo 
27T 27r 





10^ lO"^ 10^ 10^ 10^ 



Averaged squeezing bandwidth AcJo/(27r) (Hz) 



FIG. 6. Photon flux versus averaged squeezing bandwidth. 
Lines (i) to (iv) represent the amplitude squares used 
in the experiment. Trace (v) is the scaled photon flux of 
squeezing, <J/sq, which is calculated from Eqs. (|70p and (|72p . 
Dashed lines show the specific frequencies in the experiment, 
A, Q, Acjopo- 



Figure [6] shows (f/gq as a function of the squeezing 
bandwidth Aujq. We also plot experimental amplitude 
squares |ap = 1.02, 1.88, 2.87, 6.24 xlO^ s-^ In Fig. El 
there is a certain region which satisfies Aujq > Q^X 
and ^/gq < |<^P- For example, let's set the effective 
squeezing bandwidth as ten times of the resonant fre- 
quency, AuJq^ = lOQ (> 10 A). In this case, we obtain 
^/gq = 1-37 X 10^ s~^ which is still an order smaller than 
the experimental Thus we can assume the effec- 

tive squeezing bandwidth which simultaneously satisfies 
AcjOq^ ^}^X and ^/gq <^ |<^P- Accordingly we may con- 
clude that the approximation Sai{^) ^ |ape^^p is valid 
within our experimental conditions. 



